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Abstract. We consider the linear water-wave problem in a peri- 
odic channel which consists of infinitely many identical containers 
connected with apertures of width e. Motivated by applications 
to surface wave propagation phenomena, we study the band-gap 
structure of the continuous spectrum. We show that the for small 
apertures there exists a large number of gaps and also find asymp- 
totic formulas for the position of the gaps as e — > 0: the endpoints 
are determined within corrections of order e 3 / 2 . The width of the 
first bands is shown to be 0(e). Finally, we give a sufficient con- 
dition which guarantees that the spectral bands do not degenerate 
into eigenvalues of infinite multiplicity. 



1. Introduction 

1.1. Overview of the results. During the past decades the physical 
aspects of wave propagation in periodic media has been extensively 
investigated because of applications of these structures. The range of 
applications varies from photonic crystals to meta-materials or Bragg 
gratings of surface plasmon polariton waveguides, from energy harvest- 
ing (piezoelectric materials) to surface wave propagation in periodic 
channels as in this paper. Although the various applications are physi- 
cally distinct, the mathematical representation of these models in their 
different contexts has many similarities. 
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Early on it was noticed that waves propagating in periodic media 
have spectra with allowed bands separated by forbidden frequency gaps. 
This phenomenon was first discussed by Lord Rayleigh [27] . It has also 
attracted some interest in coastal engineering because it provides a pos- 
sible means of protection against wave damages [151 HB], f° r example 
by varying the bottom topography by periodic arrangement of sand- 
bars. The existence of forbidden frequencies is conventionally related 
to Bragg reflection of water waves by periodic structures. Bragg re- 
flection in this context is an enhanced reflection which occurs when 
the wavelength of an incident surface wave is approximately twice the 
wavelength of the periodic structure. For this mechanism to be ef- 
fective, the waves must be relatively long in order that they can be 
influenced by changes in depth [T8] . 

A similar occurrence may also happen when waves are propagating 
along a channel with periodic variations in the channel width. In [13], 
and later [IT], the authors studied a channel, which had a wall with 
a periodic stepped structure. By the resonant interaction theory they 
were able to verify that significant wave reflection could occur. The 
above discussions are all based on the assumption of small wall irregu- 
larities. 

On the other hand, gaps in the continuous spectrum for equations or 
systems in unbounded waveguides have been studied in many papers, 
for example in [3], [I], [6], [8], [9], [21], [22] to mention some. In [25] the 
authors studied the linear elasticity system and proved the existence 
of arbitrarily (though still finitely) many gaps, the number of them 
depending on a small geometric parameter; this result is analogous to 
Corollary 13.21 below. In the setting of the linear water-wave problem , 
spectral gaps have been studied in [T2] , [TH] , [23] and [7J , though from 
a different point of view than in the present work. 

In this paper we consider the surface wave propagation via the linear 
water wave equation with spectral Steklov boundary condition on the 
free water surface, see the equations fll.8l) - fll.lOI) . which are called the 
original problem here. The water-filled domain IF forms an unbounded 
periodic channel consisting of infinitely many identical bounded con- 
tainers connected by apertures of width e > 0, see Figure 11.11 The 
first results, Theorem 13.11 and Corollary 13.21 show that the continuous 
spectrum a of the original problem (which is expected to be non-empty 
due to the unboundedness of the domain) has gaps, and the number of 
them can be made arbitrarily large depending on the parameter e. An 
explanation of this phenomenon can be outlined rather simply using the 
Floquet-Bloch theory, though a lot of technicalities will eventually be 
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involved. Namely, if e = 0, the domain becomes a disjoint union of in- 
finitely many bounded containers, and the water-wave problem reduces 
to a problem on a bounded domain (we call it the limit problem) , hence 
it has a discrete spectrum consisting of an increasing sequence of eigen- 
values (A^)^. On the other hand, for e > 0, one can use the Gelfand 
transform to render the original problem into another bounded domain 
problem depending on the additional parameter 77 G [0,27r). For each 
fixed 77 this problem again has a sequence of eigenvalues (A^.(^))^ =1 . 
Moreover, by results of [TS], [23], Theorem 3.4.6, and [20], Theorem 
2.1, the essential spectrum a of the problem f ll .8j) ( JTTTU]) equals 

00 

(1.1) ^=|J T ^ T * = { A ^) : ^[0,2vr)}, 

k=l 

where the sets T e k are subintervals of the positive real axis, or bands of 
the spectrum. (For the use of this so called Bloch spectrum in other 
problems, see for example [8], or [2].) In general, those bands may 
overlap making a connected, but in Theorem 13. II we obtain asymptotic 
estimates for the lower and upper endpoints of T|: we show that A° < 
A|(?7) < A° + C k e for all k and rj and for some constants C k > 0. In 
view of (II. ip this implies the existence of a spectral gap between T\ 
and Tk+i for small e and k such that A° 7^ A° +1 . However, since the 
estimates depend also on k, we can only open a gap for finitely many 
k, though the number of gaps tends to infinity as e — > 0. 

The asymptotic position (as e — > 0) of the gaps is determined more 
accurately in Theorems 13.51 and 13.61 those main results state that 

T £ fc = (A° fc + A k e + 0(e 3 / 2 ), A° + B k e + 0(e 3 / 2 )) 

where the numbers A k < B k depend linearly on the three dimensional 
capacity of the set 6. This result also ensures that in case A k 7^ B k 
the bands T^, do not degenerate into single points, which means that 
the spectrum of the original problem indeed has a genuine band-gap 
structure. Facts concerning the numbers A k , B k are discussed after 
Theorem 13.61 

As for the structure of this paper, we recall in Section 11.21 the exact 
formulation of the linear water-wave problem, its variational formu- 
lation as well as the parameter dependent problem arising from the 
Gelfand transform, and the limit problem. Section [2] contains the for- 
mal asymptotic analysis which relates the spectral properties of the 
original problem with the limit problem and which is rigorously justi- 
fied in Secion [3j The main results, Theorems 13. 1[ 13.51 and 13.61 as well as 
Corollary 13.21 are also given in Section [3j The proofs are based on the 
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Figure 1.1. Side view of the waveguide (a) and of the 
periodicity cell (b) 

max-min principle and construction of suitable test functions adjusted 
to the geometric characteristics of the domains under study. 

Acknowledgement. The authors want to thank Prof. Sergey A. 
Nazarov for many discussions on the topic of this work. 

1.2. Formulation of the problem, operator theoretic tools. Let 

us proceed with the exact formulation of the problem. We consider an 
infinite periodic channel IP (see (I1.7P ). consisting of water containers 
connected by small apertures of diameter 0(e). The coordinates of the 
points in the channel are denoted by x = (xi,X2,x 3 ) = (y\,y2,z) = 
(y,z), and x' = (x 2 ,x 3 ) stands for the projection of x to the plane 
{xi = 0}. We choose the coordinate system in such a way that the 
axis of the channel is in xi-direction and the free surface is in the plane 
{x 3 = 0}. 

Definition 1.1. We describe the geometric assumptions on the pe- 
riodicity cell in detail, as well as some related technical tools includ- 
ing the cut-off funcions. Let us denote by w, C M 3 a domain with 
a Lipschitz boundary and compact closure such that its intersections 
with {xi = 0}- and {x\ = l}-planes are simply connected planar do- 
mains with positive area and contain the points P° = (0, P2, P3) and 
P 1 = (1, P 2 , P3) with P3 < 0, respectively. Then the periodicity cell 
and its translates are defined by setting (see Figure II. ip 
(1.2) 

w = {x G zu, : x 3 < 0, Xi G (0, 1)}, zuj = {x : (xi—j, x 2 , x 3 ) G to}, j G Z. 

Furthermore, we assume that the set C I 2 is a bounded planar 
domain with smooth boundary 89, containing the origin. We assume 
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that 9 is so small that the set {0} x {29 + P 3 ) is contained in dzu and 
sup (x 3 + P 3 ) =: d e < 0. We define the apertures between the 

(x2,X3)e0 

container walls as the sets 

(1.3) 6) = {x = (j,x') : e" 1 ^' - (P 2 ,P 3 )) e 9}, J G Z. 

It is plain that rr 3 < for x G #j for all < e < 1, by the choice of dg. 
We will need at several places a cut-off function 

(1.4) X e G C °°(M 3 ), 

which is equal to one in a neighborhood of the set {0} x 6 and vanishes 
outside another compact neighborhood of {0} x 9. More precisely, we 
require that 

(supp ( X e) + (0, 0, P 3 )) n { Xl = 0} C dw, 

(1.5) (supp (xe) + (0, 0, P 3 )) H {xi > 0} C w 

(this is possible by the specifications made on 9) and xe vanishes, if 
> 1/4 or a; 3 + P 3 > c? e /2. We also assume that d xi xe = 0, when 
X\ = 0. Furthermore, denoting Xj( x ) = Xe{ x ~ P 5 )) it follows from 
the above specifications that Xj^) = 0, if x 3 > o?e/2; in particular Xj 
vanishes on the free water surface 7. Finally, we will need the scaled 
cut-off functions 

(1.6) X j = xe{e-\x-P i )). 

It is plain that also Xj vanishes on 7 for < e < 1 and that Xj(x) = 1 
for x e 9 e j, j = 0, 1. 

Definition 1.2. The periodic water channel is defined by 

(1.7) ir = U(o7 i u^), 

and it will be the main object of our investigation. The free surface 
of the channel is denoted by T e = <9LT D {x 3 = 0}, and the wall and 
bottom part of the boundary is S e = <9IP \ T £ . The boundary of the 
isolated container zu, the periodicity cell, consists of the free surface 7 
and the wall and bottom <j e with two apertures 9q and 9\. 

Remark 1.3. We will use the following general notation. Given a 
domain S, the symbol (•, •)= stands for the natural scalar product in 
L 2 (S), and H k (E), k e N, for the standard Sobolev space of order k 
on S. The norm of a function / belonging to a Banach function space 
X is denoted by ||/;X||. For r > and a e R N , B r (a) (respectively, 
S r (a) ) stands for the Euclidean ball (resp. ball surface) with center 
a and radius r. By C,c (respectively, C^, c^, C(k) etc.) we mean 
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positive constans (resp. constants depending on a parameter k) which 
do not depend on functions or variables appearing in the inequalities, 
but which may still vary from place to place. The gradient and Laplace 
operators V and A act in variable x, unless otherwise indicated. 

In the framework of the linear water-wave theory we consider the 
spectral Steklov problem in the channel IP, 

(1.8) -Au e (x) = 0, xGlT, 

(1.9) d n u e (x) = 0, x E S e , 

(1.10) d z u e (x) = AV(x), x G P, 

u e is the velocity potential, \ e = g~ x u 2 is a spectral parameter with u > 
as the frequency of harmonic oscillations, and g is the acceleration 
due to gravity. By the smoothness assumptions, the normal derivative 
d n is defined everywhere on <9IT except at the boundary edges. It 
coincides with d z = d/d z at the free surface Y e . 

The rest of this section is devoted to presenting the operator theoretic 
tools which will be needed later to prove our results: Gelfand transform, 
variational formulation of the boundary value problems, and max-min- 
formulas for eigenvalues. The spectral problem (jl.8j) — (jl.lOj) can be 
transformed into a family of spectral problems in the periodicity cell 
using the Gelfand transform. We briefly recall its definition: 

(1.11) v(y,z)^ V(y,z,T]) = -J=^2exp(-iri(z + j))v(y,z + j), 



2tt . 

where (y, z) G IT on the left while r\ G [0, 2tt) and (y, z) G w on the 
right. As is well known, the Gelfand transform establishes an isometric 
isomorphism between the Lebesgue spaces, 

L 2 (IT) ~ L 2 (0,2vr;L 2 (ro)), 

where L 2 (0,27r; B) is the Lebesgue space of functions with values in 
the Banach space B endowed with the norm 

\\V-L\^-B)\\= n 27r \\V(r)y,B\\ 2 dri\ ' . 

The Gelfand transform is also an isomorphism from the Sobolev space 
H l (Il e ) onto L 2 (0, 2vr; H l e>n {w)) for I = 1, 2. The space H^w) consists 
of Sobolev functions u which satisfy the quasi-periodicity conditions 

(1.12) w(0, x') = e- ir >u(l, x'), (0, x') G 9q, 

(1.13) d Xl u(0,x') = e-^d Xl u(hx') } (0,x) G 9 e Q , 

whereas H\ Jvj) is the Sobolev space with the condition (I1.12p only. 
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Applying the Gelfand transform to the differential equation ( 11. 8ft and 
to the boundary conditions (I1.9I) - (I1.10I) . we obtain a family of model 
problems in the periodicity cell w parametrized by the dual variable 77, 

(1.14) - AU e (x;r]) = 0, xEzu, 

(1.15) d n U e (x;r]) = 0, xEa e , 

(1.16) d z U e (x;rj) = A e (r/)[/ e (x; r/), x G 7, 

(1.17) f/ e (0, a;'; 77) = e~ ir, U € {l, x'\ 77), iG 

(1.18) £^(0, 2'; 77) = e-^tf^l, x'; 77), xG^. 

Here, A 6 = A e (?7) is a new notation for the spectral parameter X e . 

The apertures disappear at e = so in that case the also quasi- 
periodicity conditions cease to exist. Hence, we can consider the prob- 
lem ( ll.14p - fl.18p as a singular perturbation of the limit spectral prob- 
lem 

(1.19) -AU°(x) = 0, xEzu, 

(1.20) d n U°(x) = 0, x E a, 

(1.21) d z U°{x) = A°U°{x), 2G7 

with A clS db spectral parameter. 

The variational form of the problem ( ll.14p - fl.18p for the unknown 
function U e E H\ (w) can be written as 

(1.22) (Vf/ e , W) ro = A e (U% V)^VE Hl v (w), 

and the corresponding variational formulation of the limit problem for 
U E H 1 (vj) reads as 

(1.23) (W, VV% = A(U, V) 1 , V E H\w). 

We denote by T-L e the space H\ Jfjj) endowed with the new scalar 
product 

(1.24) (u,v) £ = (Vu,Vv) m + (u,v) y , 

and define a self-adjoint, positive and compact operator £> e (r/) : H e — > 
T-L e using 

(1.25) (£ £ (r/) M ,t;) e = (u, W ) 7 . 

The problem (11.221) is then equivalent to the standard spectral problem 

(1.26) B e (r])u = M e u 
with another spectral parameter 

(1.27) Af e = (1 + A e ) _1 . 
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Clearly, the spectrum of B e (r]) consist of and a decreasing sequence 
(M^(r/))^ =1 of eigenvalues, which moreover can be calculated from the 
usual min-max formula 

(1.28) M e k ( V ) = minmax^^, 

where Ef. C "H e is any subspace of co-dimension k — 1. Using (ll.24p 
and fll.25p . we can write a max-min formula for the eigenvalues of the 
problem (TE22jl : 

* e / \ 1 , ■ (Vv, Vw) ro + (v,v) 7 

^k\V) = ~~ = maX 111111 



(1.29) = max mm — — 

£ fc veE k \\v; L 2 (7)|| 2 

On the other hand, the connection ( 11.271) and the properties of the 
sequence 

{M^(r])) k % 1 mean that the eigenvalues (11.291) form an unbounded se- 
quence 

(1.30) < A\(r,) < Al(r,) <...< Al( V ) < . . . -+ +oo. 

The eigenfunctions can be assumed to form an orthonormal basis in the 
space L 2 {w). The functions r\ i— > A.%(r)) are continuous and 27r-periodic 
(see for example [TO] . Ch. 9). Hence the sets 

(1.31) r h = {Al(r,):r,€[0,2n)} 

are closed connected segments, which may degenerate into single points; 
their relation to the original problem was already mentioned in (II. ip . 

For the same reasons as above, and working in the space H x {vj) in- 
stead of H e (with norm induced by (jl.24j) . which is equivalent to the 
original Sobolev norm of H 1 (vj)), the limit problem with e = has 
an eigenvalue sequence (A°)^ =1 like (11.301) . however, the eigenvalues 
are constants with respect to the ^-parameter. The first eigenvalue 
equals 0, and the first eigenfunction is the constant function. Analo- 
gously to (11.291) we can write 

(1.32) A° = max min ^'^/T^ > 

F k veF k ||t;;L 2 (7)|| 2 

where again Fk C i7 1 (ro) is an arbitrary subspace of codimension k — 1. 
We denote by 

(i-33) (f4°)r=i 

an L 2 (7)-orthonormal sequence of eigenfunctions corresponding to the 
eigenvalues ( 11.321) . 
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Lemma 1.4. For all k there exists a constant C k > such that 

(1.34) \U° k (x)\<C k , \VU° k (x)\<C k 

for all x G supp (xj)r\m, j — 0, 1 (and hence for all x G supp (Xj)r\w, 
<e<l). 

Proof. Let for example j = (the other case is treated similarly), 
and define the domains G\,G2 C M 3 with C°° boundary such that 
Go : = supp (xj) C G\ C Gi C C7 2 C {x 3 < 0} and G2 still so small 
that 

(1.35) G 2 H {x x = 0} C dw and G^n {xi > 0} C w. 

As a consequence, these domains are smooth enough so that we can use 
the local elliptic estimates [T], Theorem 15.2, to the solutions U® of the 
equation (II. 19ft : this yields for every / = 1,2, . . ., a constant > 
such that 

\\U° k -H l+1 (G n nzu)\\ <c l 4\\u!-,H l -\G n+1 nuj)\\ + \\u!-L\G n+1 nuj)^ 

for n = 0, 1. Applying this first with n = 1 and I = 1 we get a 
bound for \\U k ; H 2 (G\ fl zu)\\ and then, with n = and Z = 2, for 
||^;if 3 (G nci7)||. The standard embeddingsiJ 2 (Gino7) C C B (G Q nw) 
and iJ 3 (G n w) C G^(G n ro) imply the result. □ 

2. The formal asymptotic procedure 

2.1. The case of a simple eigenvalue. To describe the asymptotic 
behavior (as e — » 0) of the eigenvalues A|(^) of the problem (11.141) - 
(ll.lSp we consider first the case A° is a simple eigenvalue of the problem 
(jl.l9p ~ (jl.2ip for some fixed k. Let us make the following ansatz: 

(2.1) At( V )=A° k + eA' k ( V )+Al( V ), 

where A' k {ri) is a correction term and A k {rj) a small remainder to be 
evaluated and estimated. While the remainder is treated in Section 
13.21 we derive here the expression (I2.13P for A' k (rj), cf. also (12.181) and 
dZA9]). 

The corresponding asymptotic ansatz for the eigenfunction reads as 
follows: 

(2.2) W k (x; V ) = U° k (x) 

+ Xo(x)w ko (e- 1 (x - P )) + X i(x)w kl (e^(x - P 1 )) 

+ eU' k (x;r]) + U e k (x; rj), 

where (U k ) k L 1 is as in (jl.33p . The functions w k Q and w k i are of bound- 
ary layer type, and Xj is given above (11.61) . 
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The boundary layers w k j depend on the "fast" variables ( "stretched" 
coordinates) 

e = (&,&m=e-\x-Pi), j = 0,1. 

They are needed to compensate the fact that the main regular term U k 
in the expansion does not fulfil the quasi-periodicity conditions ( I1.12p - 
( I1.13P . By Lemma 11.41 and the mean value theorem, the eigenfunction 
U®(x) has the representation 

U° k (x) = U° k (n + 0(e), xe6] 

near the points P\ j = 0, 1. Hence, w k o and w k \ should be the solutions 
of the problems 

A 5 o Wfeo (£°) = 0, f? > 0, 
%^ (O = 0, ee{0}x(R 2 \9), 

w k0 (e) = a fe o, £° e {0} x 9, 

and 

A e w kl (e) = 0, £<0, 

d e w kl (e) = o, e 1 e{o}x(M 2 \^), 

w kl (^) = a kl , G {0} x 9 

in the half spaces > 0} and {£\ < 0}, respectively; the numbers 
Qjfzj are explained below. Both of the functions w kj , j — 0,1, can be ex- 
tended to even harmonic functions in the exterior of the penny-shaped 
domain {0} x 9: 

(2.3) A. w kj (?) = 0, f e K 3 \ ({0} x 6), 

w kj (e) = a kj , ? ed({0}x6). 
Furthermore, the problem (12.31) admits a solution (see 

cap 3 9 

w 

(2.5) w k 0) = o(ier 2 ) , v £! /a,(e) = o(ier 3 ), 

where cap 3 (#) is the 3-dimensional capacity of the domain {0} x 9 and 
(12.51) concerns large e-behavior. Moreover, the solution has a finite 
Dirichlet integral: 

(2.6) / |v^(e)| 2 de<o 

for some constant C > 0. 



(2.4) Whj ^) = ak .-^- + yj k .^ 
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The coefficients a k j can be found from the quasi-periodicity condi- 
tions ffTTITjl - ffTTTgj) . Clearly, for each e > 

Ul(P°;r ] ) = e-^ k (P 1 ; V ), 

which together with the asymptotic expansions yield the relations 

U° k (P°) + a k0 = e^iU^P 1 ) + a kl ) and a k0 = -e~^a kl 
for the coefficients. Hence, 

(2.7) akl = -e^a k0 , a k0 = ± (e^U^P 1 ) - U° k (P )) . 

Now we can write a model problem for the main asymptotic correc- 
tion term U k . 

(2.8) - AU' k (x; rj) = AW k (x) xEzu, 

(2.9) (d z - A° k )U' k (x; rj) = K(v)U° k (x), x e 7 , 

(2.10) d n U' k (x;rj) = 0, xGct, 
where we denote 

(2.11) w k (x) = ^xj(x) ^^ y xe*,,. 

The number A.' k (j]) will be unique for given k and rj; this fact follows 
from the solvability condition in the Fredholm alternative, see Lemma 
12.11 and its proof below. It must be determined as follows by using 
the Green formula and the normalization in fll.33[) ( ds is the surface 
measure) : 

K(v) = K(v)\\U k ;L 2 ( 1 )f = 



(d z U' k (x; rj) ~ A° k U' k (x; rj)) U° k (x) ds{x) 



U° k (x)d n U' k (x; rj) - U' k (x; V )d n U° k (x) ds(x 



(2.12) = / U°(x)AU' k (x]rj) = - / U%(x)AW k (x)dx. 



Taking into account that the last integral converges absolutely and 
using the Green formula again yield 

K(V) = limE / Wx)9 n (- a -f^) ds(x) 

r ^°j^Jsr{Pi)n™ V \x-P3\ J 



27rcap 3 (e) [a k0 U°(P°) + a kl U^{P r 
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see (12. lip and Remark 11.31 for notation. According to ( 12.7ft we finally 
obtain 

(2.13) A' k ( V ) = nc^ 3 (6)\U° k (P ) - e-^P 1 )! 2 . 



Lemma 2.1. Choosing A' k (rj) as in (12TT3|) . the problem (12T8|) - d2TT0|) 

has a solution U' k G H 1 ^). 

Proof. The variational formulation of the problem (I2.8p — (I2.10p reads 
as 

(2.14) (VC/£, W) ro - A° k (U' k , V), = (VWik, W) ro - A' fc (I& V) 7 . 

We remark that the function l/\x — Pi\ is harmonic in w, and since 
Xj equals constant one in a neighborhood of P\ the function AW k 
vanishes there, hence, W k and VW k are smooth as well as uniformly 
bounded everywhere in w. Moreover, U k G L 2 {j). 

Using the definition fl 1 . 2 5 [) of the operator B = B°(rj) which does 
not depend on 77 and acts Sobolev space H l {w) (with norm induced 
by (11.241) ) we can rewrite (12. 14ft as follows: 
(2.15) 

(U' k , V) - (A° k + l)(BU' k , V) = (W k , V) - A' k {BUl V) - (BW k , V) 

which means that U k must be a solution of equation 

(2.16) (B - M° k )U' k = -M° k (W k - BW k - A' k BU° k ). 

Notice that U k is the solution of homogeneous problem (12.161) . so, by 
the Fredholm alternative, (I2.16P is solvable, if and only if the right hand 
side of it is orthogonal to the function U k . This condition is satisfied 
by choosing A k (rj) as above, since 

(W k - BW k - A' k BU° k , U° k ) = (W*, Vf/ fc °) ro - A' k \\U° k , L 2 ( 7 )f = , 

by (12TT2D and (VW fc , Vf/ fe °) ro = -(AW k ,U^)„. This last identity fol- 
lows from the first Green formula, because the normal derivative of 
W k vanishes on dw due to the properties of the function Xji see below 

(USD. ' ' □ 

2.2. The case of a multiple eigenvalue. In this section we complete 
the asymptotic analysis by studying the behavior of eigenvalues AK77) 
in the case some A k has multiplicity m greater than one: we have 

A < A - - A < A 

The ansatz (12. ip is used again. Furthermore, in accordance with (11.331) 
we denote by {U k+ j)o<j< m -i C L 2 {^/) an orthonormal system of eigen- 
functions associated with the eigenvalue A k . Any eigenfunction U° 
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corresponding to A° can be presented as a linear combination 

171 — 1 

3=0 

Analogously to (12. 2p we introduce the asymptotic ansatz 

(2.l7p e (x;ri) = U°(x) 

+ Xo{x)w ko {e-\x - P )) + X i(x)w kl (e-\x - P 1 )) 
+ eU'(x;rj) + U e {x;r]). 

Using the same argumentation as in the previous section we construct 
the boundary layers w k j, j — 0,1, which satisfy the conditions 

^(e) = a kJ c -^ + o(\e\- 2 y, 

here the coefficients a k j come from the equations ( 12. 7ft . where U® is 
replaced by U°. The main asymptotic term U' is also treated in the 
same way as in Section [2TTI To use the Fredholm alternative for finding 
K+j(v), j = 0, • • • , m - 1, we write 

and making use of the Green formula as above we get 

m— 1 



where 



Pij = vrcap 3 (£)(f4 +i OP°) - e^UUP 1 ))^^ ) - ^^(P 1 )). 

Hence, A' k+ j(rj) is an eigenvalue of the matrix B = (/9y)r}=o- This ma- 
trix has rank one, because it can be represented in the form B = vv T , 
where v is a vector with components Vj + \ = U® + j(P°) — e~ ir, U k+ AP x ) , 
j = 0, . . . , m — 1. This means that 



m— 1 

(2.18) A' k ( V ) = 7rcap 3 (0) £ |^(P°) - e^L/j^ 1 )! 2 , 



(2.19) A' (7/) = 0, l<j<m-l 
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3. Justification of asymptotics 

3.1. Existence of gaps. The first estimate on the eigenvalues of the 
problem fjl.22j) can now be stated as follows. 

Theorem 3.1. For any k G N there are numbers tk > and Ck > 
such that for every e G (0, e^) and any dual variable rj G [0, 27r), 
the eigenvalues of the problem (11.221) and the eigenvalues of the limit 
problem (I1.23P are related as follows: 

(3.1) A° < A%( v ) < A° k + C k e. 

As mentioned in the introduction (see the explanations around (11.11) 
and (11.311) ). this result implies the existence of any prescribed number 
of gaps in the continuous spectrum a, since (13. ip also establishes an 
estimate for the endpoints of the intervals T^,: one needs to take enough 
many distinct eigenvalues A° and a small enough e. 

Corollary 3.2. Given any number N G N, the continuous spectrum 
a of the problem (I1.8l) - (ll.l0p on IP has at least N gaps, if e is small 
enough. 

Proof of Theorem \3.1[ We apply the max-min-principle described in 
Section 11.21 and first prove the estimate 

(3.2) AJfo) > A°. 

Indeed, we recall that in the equations (11.321) and (ll.29p both Ff. C 
H x {w) and Ek cK e = H^ v (zu) are arbitrary subspaces of co-dimension 
k — 1. Since H^^zu) C H l {w), each Ek is contained in some Fk, and 
thus the infimum in (ll.32p is smaller than that in (ll.29p . 

So we turn to the upper estimate in (13.11) and fix a k G N. Let the 
eigenf unctions U® be as in (I1.33P and let Hk C H e (w) be a subspace 
spanned by the functions Y e Uj, where j = 1, . . . , k and 

(3.3) Y e = 1 -X -X x G C°°(w) 

and Xj are as in (jl.6p . We remark that the functions Y £ Uj satisfy the 
quasiperiodicity condition fll. 12[) in the definition of the space H e , since 
Y e vanishes in a neighborhood of the apertures, see the remarks around 
(11.61) . Moreover, the sequence (YJJ® , YJJ®, ■ ■ ■ , Y t U®) is still linearly 
independent, due to the L 2 (7)-orthogonality in (I1.33P and the fact that 
Y e equals 1 in the set 7. Hence, the dimension of Hk is k. 

If Ek is an arbitrary subspace of H e of co-dimension k — 1 (cf. (I1.29P ). 
the intersection Ek PI Hk contains a non-trivial linear combination 

k k 

(3.4) U(x) = Y e (x)J2 a i U j( x )> Yl l^l 2 = L 

3=1 j=i 
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By the remarks just above we have \\U; £ 2 (7)|| = 1. Hence, from fll.29|) 
and (I1.33P we infer that 



|W; L 2 (zu) 

\\U;V{-yW 



\VU: Lriw) 



+ 



K K 

3=1 3=1 

k 

V(^a j f/°);L 2 M|| 2 + 2(V((-X -X 1 )^a,^) ) v(^a^° 

3=1 

k 

v(3(JCb + X 1 )^a^);L 2 ( ro ) 



3=1 



3=1 



3 = 1 



To evaluate the first term on the right hand side notice that the func- 
tions Uj satisfy (I1.23P so that the L 2 (7)-orthogonality of (11 .33j) implies 



3=1 



3,1=1 



(3.6) = £ C/°) 7 = s 2A ? < A °> 

3,1=1 3=1 

where the last inequality follows from H3.4[) and the fact that the eigen- 
values A° are indexed in increasing order. Furthermore, we use Lemma 
11.41 as well as the facts that the supports of Xi, I = 0, 1, have measure 
of order e 3 , IVX/I are of order e _1 , and \aA < 1 to estimate 



V ((X + Xx) £ *jUj) . V ( Yl a i U j 

3=1 3=1 

< k 2 (sup(l,\U°(x)\,\VUHx)\)) 2 sup(l,|VX(x)|) / dx<C k e 2 , 

h 2 

\v( K (x +x 1 )j2a 3 u!;y,L 2 (w)y 

3=1 

(3g|c 2 (sup (1, |Vf/°(x)|)) 2 (sup (1, |VX(x)|)) 2 j dx < C k e, 



where S = supp(X + Xi). Combining this with (13 ,6p and (13. 5 p yields 
the result. □ 
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3.2. Asymptotic position of spectral bands. In this section we 
will prove the validity of the asymptotic ansatz (12 .ip , see Theorem 13.51 
This yields our main result concerning the asymptotic position of the 
spectral bands, Theorem 13. 61 

We start the proof by recalling a classical lemma on near eigenvalues 
and eigenvectors (see [28] and also, e.g., [5]). 

Lemma 3.3. Let B be a selfadjoint, positive, and compact operator 
in a Hilbert space %. If a number fi > and an element V G 7-L 
satisfy ||V;H|| = 1 and \\BV — fiV;7i\\ = r G (0,fi), then the segment 
\p — r, fi + r] contains at least one eigenvalue of B. 

To apply Lemma [3.31 to the operator B e {ri) of (ll.25p . we fix an arbi- 
trary k and, keeping in mind the formula (11.271) . define the approximate 
k:th eigenvalue and eigenvector of B e (rj) by 



(3.8) 

where 
(3.9) 



fi k = (1 + A° + eA'^))- 1 , 
V fc (x) = WUu-WW^Ukix), 

U k (x) = (l-X (x)-X 1 (x))U%(x) 
+ X (x)U° k (P°) + X^U^P 1 ) 
+ Xo{x)w kQ {t-\x - P )) + xi{x)w kl {e-\x - P 1 )) 
+ e(l-X (x)-X 1 (x))U' k (x,r l ), 



Ul is as in (OB]) . Xj(x) = X e(e 



[x 



P- 7 )) and Xj are as i n (H-61) . 



We need a lower bound for the norm of XA k . 
Lemma 3.4. For all k there exists a constant C k > such that 



(3.10) 



e II 2 



< C k e^ . 



Proof. Recall that the expression for U k , (I3.9p . contains the term U k ; 
let us denote U k := U k - U%. By Ol . f02lH03D . (Q5]l . we have 
||[/°;ft e || 2 = ||f/ fc ;L 2 ( 7 )|| 2 + ||W°;L 2 (w)|| 2 = 1 + A°. Hence, by the 
Cauchy-Schwartz inequality, 

\\u k -w\\ 2 -\-^ k 



2(U° k ,U k ) e + \\U k ;n 



e II 2 



(3.11) < 2^/1 + A°||Z4; W\\ + ||Z4; ft £ || 2 . 

Taking into account the definition of the norm of H 6 , the formula 
(13. 9 p and the fact that the functions Xj an d Xj vanish on 7 we find 
that \\U k ] "H e || is bounded by the sum of the expressions 



(3.12) 



\v(xM-u°k{P j )))-,L 2 {™) 



0,1, 
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(3.13) WVixjWkAe'Hx-P^^WW J = 0,1, 

(3.14) \\e(l-X -X 1 )U k ;H 1 (w)\\. 

First we use the observation (Lemma 11.4ft that the supports of the 
functions Xj and VXj are contained in balls of radius 0(e) and that 
|f/°(x)| and |Vt/°(x)| are uniformly bounded in these balls, hence 
U k (x) — U k (Pi) = 0(e) there. So, (13.121) can be bounded by a constant 
times 

/ \U° k -U° k (P j )\ 2 \VX j \ 2 dx + [ \VU° k \ 2 \X 3 \ 2 dx 

J zu J zu 

(3.15) <C([ e 2 e~ 2 dx + f dx)^ < Ce 3 ' 2 . 

We estimate the terms (13. 13ft using the fact that the support of the 
function Vxj is contained in a set {c < \x — P J '| < C} =: Sj for 
some constants < c < C (see above (11.61) ). hence, by the estimate 

(3.16) |wjy(e _1 (a: - P j ))\ < Ce\x - P J \- 1 for x E Sj. 
Applying (12.61) yields 

\V{x 3 {x)w k] {t-\x-P' J )))\ 2 dx 
I \w k3 (e~\x-PJ))\ 2 dx+ I \V{w k] {e-\x-P>)))\ 2 dx 

J Sj J ZU 

< [ Ce 2 \x- P j \- 2 dx + e\- 2 [ \Vzjw kj (Z j )\ 2 dt j 

< C x e. 

Finally, by Lemma [2.11 U' k belongs to the space H 1 (w). For the terms 
(I3.14p we thus get the bound 

\\e{l-X -X 1 )U' k ;H 1 (w)\\ 

< Ce\\U' k -H 1 (w)\\ +Oemax||VX,-;L 2 (u7)|| 1/2 ||^;L 2 (u7)|| 1/2 < O'e, 

i=o,i 

since || VX,-; L 2 (ot)|| < Oe 1//2 , due to the measure of the support of 
VXj. □ 

As a corollary of this lemma, if e 6 (0, eo], then the bounds 

(3.17) 0</i fc <c M \\U;H e \\ > c u > 0, 

hold true with some positive constants c M and cu depending on w and 
9 only. 



< 
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The next theorem provides quite accurate asymptotic information 
on the eigenvalues of the model problem and in particular justifies the 
ansatz (12. ip . 

Theorem 3.5. For every k > 1 there exists a constant C k such that 
for each 77 G [0, 2n) 

(3.18) |AJfo) - A£ - eA'Ml < C fe e 3/2 • 

where A' k (r]) = it cap 3 (6)\U k (P°) — e~ ir? f/°(P 1 )| 2 in the case of a simple 
eigenvalue (cf. (12 . 1 3[) ) and A' k {rj) is given by the formulas (I2.18p - (I2.19I) 
in the case of a multiple eigenvalue. 

Proof. We apply Lemma [3.31 to the operator B e (r]) with [i and V as in 
(13. 8p . Our aim is to show that r of the lemma can be chosen as small 
as C^e 3 / 2 The lemma then gives an eigenvalue M(e,n) of B e {rf) with 
the estimate 

\M(e,r)) - v k \ <C k e 3 / 2 . 

Using (13 .8p and (I1.27P this turns into an eigenvalue A(e, n) (of (I1.14p - 
(I1.18P ) satisfying (I3.18P in the place of A^ry) . However, if e is small 
enough, Theorem 13.11 guarantees that in a neighborhood of A k there 
is only one eigenvalue of the model problem, namely A\(rj). So A(e,7/) 
must coincide with it, and the estimate (I3.18P follows. 
We are thus left with the task of proving 

r=\\B\r,)V k -^ k V k ;W\\ <C k e 3 ' 2 . 

To this end we write, using V k = c^U k , (TL25|) . (Olj) . fl3T8|) . flSTTTj) . 

r = sup \{B e {r])V k - /i k V k , Z) e \ 
z 

= c u ] sup \ (U k , Z) n - fi k (U k , Z) 7 - /ifc(VW fc , VZ) J 

z 

(M%c^ sup I (A°. + eA' k )(U k , Z\ - (VZ4, VZ)J =: c^ 1 sup \T(Z)\ . 

z z 

The supremum is calculated here over all functions Z e "H e with unit 
norm. The expression T(Z) can be represented as a sum of the terms 

S 1 (Z) = -(Vf/°,VZ) ro + A°(f4°,Z) 7 , 

5 2i (z) = - (v(x,-(t/£(i»") - £4°)), vz) ro , 

S 3 (Z) = -( V (xo^o(e- 1 (x-P°)) + OT i(e- 1 (a;-P ))),VZ) ro 

+ e(Vu> fc ,VZ) ro , 
5 4 (Z) = -e(yU' k ,VZ) m -e(Vw k ,VZ) a 

+ eA' k (UlZ), + eA° k (U' k ,Z), } 
S B (Z) = -e(V((A + X 1 )U' k ), VZ) ro + e 2 A' k (U' kl Z) T 
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where w k is given by 

(3.20) ai{x) = Xo{x fJ^^ + Xl{x fJ^>M. 

First we note that S\(Z) = S±(Z) = 0, because U® and U' k are the 
solutions of the problems fO^ - fOT]) and (jX5]l -( E30]l . respectively; 
see also (gZHJ) . 

To estimate S<ij{Z) we use the Cauchy-Schwartz inequality: 

|(v(x,(^(p^)-^°)),vz) ro | 

Here the last inequality follows from an estimate already made for 
(I3.12p (see (I3.15P ) and the assumption on the norm of Z. 

The first term in S${Z) can also be treated using the Cauchy-Schwartz 
inequality, the properties of the cut-off functions Xj (once again as in 
(En3]) ) and Lemma O 

e\(V((X + X 1 )U' k ),VZ)J 

< e\\ V((X Q + X 1 )U' k ); L 2 (w)\\ \\Z; H\w)\\ 

< e ||V(X + X 1 );L 2 (u7)||||^;L 2 (u7)|| 

+ e\\X + Xi, L 2 (w)\\\\VU k ; L 2 (m)\\ <C k e 3/2 . 

The surface integral in S$(Z) can be estimated simply by the trace 
inequality. 

To provide an upper bound for S 3 (Z) we notice that by (12. 4p . 

(3.21) w kj (e~\x - Pi)) = e a ^^ + w kj (e~\x - &)) , j = 0, 1, 

hence, using (I3.20p we can write 

l 

S 3 (Z) = (V( Xj w kj (e-\x - n)),VZ)^. 

3=0 

After integrating by parts and taking into account that w k j are har- 
monic functions we obtain 

i 

ft^ = £ faii^i? - pj )) A X3, Z) v 

3=0 

(3.22) + 2((V Xj )V(i fcj (r 1 (x-F))) 1 Z) ra 

In the second term, the support of the function Vxj is contained in 
a set {c < \x — Pi\ < C} =: Sj for some constants < c < C, hence, 
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by the estimate ( 12. 5ft . 
(3.23) 

^(wkjie^ix-P^ < Ce-^e-^x-P^y 3 = Ce 2 \x-P j \~ 3 for x G S r 
Hence, 

{(V Xj )V{w kj (e-\x-P j ))),Z)^ 
(3.24) < ( y ^(^(e-^x-P^))! 2 ^) 172 !^;^ 2 ^)!! < C'e 2 - 

The first term in ( 13. 22ft is treated with a similar argument, since <Sj 
still contains the support of Axj and the estimate 

\w kj (e-\x - P j ))\ < Ce 2 \x - P j \- 3 for x G 

again holds, by (12.51) . 
We thus get the bound 

r < C k (6)e 3 / 2 . 

□ 

As a consequence we can now provide the asymptotic widths and 
positions of the spectral bands. 

Theorem 3.6. In the case of a simple eigenvalue A k the band T e k 
of the continuous spectrum of problem (jl.8p - fll.10p has the asymptotic 
form 

= (A° fc + A k e + 0(e 3 / 2 ), A° + B k e + 0{e 3 ' 2 )) 

where 

(3.2^ = 7rcap 3 (#)min{|f4 (P°) - ^(P 1 )) 2 , \U° k (P°) + U^P 1 )] 2 }, 
(3.28,) = 7rcap 3 (e) max{|[/°(P°) - ^(P 1 )) 2 , \U° k (P°) + U^P 1 )] 2 }. 

Remark 3.7. Returning to the band-gap structure of the Bloch spec- 
trum (11.11) . the mere existence of spectral gaps does in general not 
exclude the possibility that the bands between gaps could degenerate 
into single points. This is exactly the case when the band consists of 
a single eigenvalue of infinite multiplicity. However, by Theorem 13. 6j 
if k is such that both U k (P°) and U k (P l ) are nonzero and e is small 
enough, the numbers A k and B k are distinct, and in this case the spec- 
tral band T^. does not degenerate into a point. This obviously provides 
a way to construct examples where the spectrum of the linear water- 
wave problem has a genuine band-gap structure with at least a given 
number of spectral gaps, cf. Theorem 13.21 
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